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ABSTRACT 

We have investigated the abundance anomalies of lithium for stars with plan- 
ets in the temperature range of 5600-5900 K reported by Israelian and coworkers, 
as compared to 20 normal stars at the same temperature and metallicity ranges. 
Our result indicates a higher probability of lithium depletion for stars with plan- 
ets in the main-sequence stage. It seems that stellar photospheric abundances 
of lithium in stars with planets may be somewhat affected by the presence of 
planets. Two possible mechanisms are considered to account for the lower Li 
abundances of stars with planets. One is related to the rot at ion- induced mix- 
ing due to the conservation of angular momentum by the protoplanetary disk, 
and the other is a shear instability triggered by planet migration. These results 
provide new information on stellar evolution and the lithium evolution of the 
Galaxy. 

Subject headings: planetary systems - stars: late-type - stars: abundances - 
Galaxy: evolution 



1. Introduction 



Since the first announcement of a planet orbiting a Sun-like star, 51 Peg, by Mayor 
& Queloz (1995), around 140 stars with planets (SWPs) have been reported. The main 
properties of SWPS have been investigated by many works. Among these, spectroscopic 
observations of most SWPs based on high-resolution and high signal-to-noise ratio spectra 
have been carried out. All of these works indicate a high mean metallicity of SWPs relative 
to normal field stars. Are there any other chemical signatures for SWPs? There is no 
conclusive result on this issue, and only a few works have reported abundance anomalies 
in some elements. For example, Gonzalez & Laws (2000) suggested a high [C/Fe] ratio for 
SWPs as compared with the normal stars from Gustafsson et al. (1999). However, they later 
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reported that there is no difference in abundance for C, Na, Mg, and Al elements between 
SWPs and normal stars (Gonzalez et al. 2001). Present knowledge on the abundances of 25 
elements (except for Li and Be) available in the literature shows no special feature for SWPs 
except for the high overall metallicity. 

For lithium abundance, Israelian et al. (2001) claimed that 6 Li was detected in the 
atmosphere of HD 82943, while Reddy et al. (2002) did not confirm this detection. Mandell 
& Ge (2004) reported no 6 Li presence in several lithium-poor stars. For 7 Li, King et al. (1997) 
first found several SWPs with low Li abundances, and they related them with the presence of 
a planetary companion. Gonzalez & Laws (2000) concluded that smaller Li abundances are 
found for SWPs after correcting for trends with temperature, metallicity, and age. However, 
Ryan (2000) suggested that there is no significant difference in lithium abundance for SWPs, 
and later work by Gonzalez et al. (2001) agreed with Ryan's conclusion. With a large sample, 
Israelian et al. (2004) suggested again that in the temperature range 5600-5850 K, SWPs 
show additional depletion of Li compared to normal stars. They also showed that for higher 
temperatures, the Li abundances were similar for SWPs and for normal stars. Recently, 
measurements of Be abundances by Santos et al. (2004) indicated further evidence on the 
composition difference between stars with and without planets. 

Lithium is a special element because it is easily destroyed in stars during even pre-main 
sequence and main sequence evolution, which is not the case for other elements. Therefore, 
Li abundance reflects the mixing history of the star. For SWPs, Li is further related with 
the accretion of material and the angular-momentum evolution of the system, and thus 
Li abundances have been used to distinguish the current hypotheses on planet formation 
mechanisms. In this respect, it is important to verify whether Li abundance anomalies for 
SWPs at the temperature range of 5600-5850 K reported by Israelian et al. (2004) are real, 
which will provide new information on the planet formation process and its influence on host 
stars. 

In this work we conduct a parallel abundance analysis of lithium abundances for 16 
SWPs and 20 comparison stars in the same metallicity of [Fe/H] ~ 0.1 dex at the temper- 
ature range of 5600-5900 K in order to investigate whether any anomalies persist. Further- 
more, this study is quite interesting for our understanding of the lithium behavior of solar- 
metallicity stars, which show a large scatter at 5600 < T c s < 5900 K and which the standard 
model of stellar evolution cannot explain. With a consistent analysis procedure, this work 
aims to investigate the effect of planet presence on stellar evolution. This will provide new 
information on the chemical evolution of lithium in the Galaxy. 
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2. Observations and Abundance Analysis 

SWPs are selected from a catalog 1 in the temperature range of 5600-5900 K, and 
their spectra are available from our observations with the 2.16m telescope of the National 
Astronomical Observatories (Xinglong, China) with a resolving power of 37000, or from 
Prugniel & Soubiran (2001), who observed with the 1.93m telescope of the Observatoire de 
Haute Provence with R~42 000. Prugniel & Soubiran (2001) present a database of high- 
resolution spectra of 709 stars covering a large range of atmospheric parameters, which 
enables us to select a comparison sample of stars spanning the same ranges in temperature, 
gravity and metallicity as the above sample of planet stars but not having been reported to 
harbor planets. That is, both samples cover the atmospheric parameters of 5600 < T c g < 
5900 K, log g > 3.8, and [Fe/H] > —0.3. In the selection, stromgren uvby indices from 
Olsen (1983, 1993) are adopted and the metallicity is estimated based on the calibrations by 
Schuster & Nissen (1989). Along with the initial metallicity, temperature and gravity are 
derived in the same way as described below. All spectra have signal-to-noise ratios above 
150 pixel" 1 at the Li I A6707 line. 

The spectra were reduced using standard MIDAS (for Xinglong data) and IRAF (for 
ELODIE data) routines for order definition, background correction, flat-fielding, extraction 
of echelle orders, wavelength calibration, and continuum fitting. The Fe II lines suitable for 
measurement were carefully selected, and the equivalent widths (EWs) were obtained by the 
fitting of a Gaussian function. Atomic data for Fe II lines were the same as those in Chen 
et al. (2002). The wavelengths and oscillator strengths of the Li I A6707 line are taken from 
Smith et al. (1998). 
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Table 1: Stellar Parameters for SWPs and the Comparison Stars 
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Effective temperatures were determined from the Stromgren (b — y) color index (Olsen 
1983, 1993) using the calibration of Alonso et al. (1996). The gravities were determined via 
Hipparcos parallaxes (Perryman et al. 1997), and the microturbulences were obtained by 
forcing Fe I lines with different strengths to give the same abundances. Here the oscillator 
strengths of the Fe I lines are taken from Chen et al. (2000). The new metallicity 
derived from spectroscopic analysis was updated, and the procedure of atmospheric parame- 
ter determination was iterated for consistency, even though there is no systematic difference 
between the initial metallicity and the new value, with a deviation of 0.01 ± 0.15 dex. The 
internal errors of the atmospheric parameters are around 70 K in temperature, 0.1 dex in 
log g, 0.1 dex in [Fe/H], and 0.3 kms" 1 in microturbulence. But the absolute uncertainties 
of the stellar parameters may be slightly larger. 

We have found that the photometry-based temperatures in the work are 130 K lower 
than the spectroscopically-derived values in Santos et al. (2004) for SWPs in common. 
However, our temperatures are 83 K higher than those in Takeda & Kawanomoto (2005), 
who derive temperatures by spectroscopic method for 29 stars in common. It seems that 
spectroscopically-derived temperatures have significant uncertainties, and the values depend 
on stellar model atmosphere, the selection of lines, the adopted atomic data, and the way the 
temperatures are determined. In addition, as already noted by Santos et al. (2004), surface 
gravities determined from the balance of Fe I and Fe II abundances are systematically higher 
than those derived from parallaxes (adopted in our work) by the order of 0.15-0.20 dex. 
However, iron abundances in this work and Santos et al. (2004) are similar since we adopted 
iron abundances from Fe II lines and they forced Fe I and Fe II lines to give the same 
iron abundances when they derived gravities. For the comparison stars, there are 15 stars 
in common with Edvardsson et al. (1993), and our temperatures are 87 K lower than their 
values. This is due to different temperature calibrations, and our values, based on the 
infrared flux method temperatures by Alonso et al. (1996), may be more reliable. Actually, 
the absolute temperature and gravity are not very important, since we are carrying out a 
differential analysis, and it is crucial that the stellar parameters all be in the same scale 
between the target and the comparison sample. 

The model atmospheres were interpolated from a grid of plane-parallel, LTE models 
provided by Kurucz (1995) in which convective overshoot is switched off. The ABONTEST8 
program developed by P. Magain in the Liege group was used to carry out the calculations of 
theoretical line EWs, and abundances were derived by matching the theoretical EWs to the 
observed values. Abundance errors, estimated from the uncertainties of atmospheric param- 
eters and EWs, are 0.07 dex for iron abundances, and less than 0.1 dex for Li abundances. 
Here the contribution of the 6 Li isotope to the Li I A6707 line was assumed to be negligible. 
Finally, non-LTE corrections were applied to the derived Li abundances based on the work 
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Fig. 1. — Li abundances vs. [Fe/H] for stars with (filled circles) and without (open circles) 
planets. Upper limits of Li abundances are indicated by downward-directed arrows, and the 
size of the symbol corresponds to stellar age. 
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Fig. 2. — Cumulative distribution functions based on the Kolmogorov-Smirnov test for stars 
with (solid line) and without (dash line) planets. Shown are different groups between the 
two samples, with the maximum value of D = 0.54 and the probability for the two samples 
being the same group of 0.007. 
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Fig. 3. — Stellar positions of SWPs in Israelian et al. (2004) (filled circles) with new de- 
termined temperatures as compared with Chen et al. (2001; open circles) for stars with 
5600 < T cff < 5850 K and [Fe/H] > —0.3. Downward-directed arrows correspond to stars 
with upper limits for Li abundances, and diameters of circles correspond to stellar Li abun- 
dances taken directly from Israelian et al. (2004) and Chen et al. (2001). 
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Fig. 4. — Li abundance vs. T cS . The symbols are the same as in Fig. 1. 
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Fig. 5. — Stellar positions in the HR diagram as compared with the evolution tracks from 
VandenBerg et al. (2000). The symbols are the same as in Fig. 1, but diameters of circles 
correspond to Li abundances. 
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by Carlsson et al. (1994), who studied non-LTE formation of the Li I A6708 line as a function 
of effective temperature, gravity, metallicity, and Li abundance. Stellar parameters and Li 
and Fe abundances for SWPs and the comparison stars, as well as the EWs of the lines and 
their errors are shown in Table 1. 



3. Results and Discussions 

3.1. Li Abundances versus [Fe/H] 

Li abundances as a function of metallicity for all stars in the present survey are shown 
in Fig. 1, in which stars with downward-directed arrows have an undetected lithium line, 
and the symbol size corresponds to stellar age. The upper limits of the Li abundances 
are estimated by assuming EWs of 3.0 mA, which is the estimated error in EW from the 
comparison used in Chen et al. (2000). Based on Table 1, the scatter of the EW comparison 
between this work and Takeda & Kawanomoto (2005), Chen et al. (2001), and Lambert 
et al. (1991) is around 4.0 mA, which indicates an error of 4.0/ ~ 3.0 mA assuming; the 
same uncertainties for the EWs from the other works. Note that the statistical uncertainties 
based on the formula by Cayrel (1988), presented in Table 1, are smaller. Since experimental 
errors by different authors during spectrum normalization and EW measurements are not 
considered in Cayrel's (1988) formula, we believe that an upper limit to the EWs of 3.0 mA 
for stars with an undetected Li line may be more reasonable. 

An interesting result from Fig. 1 is that there is a lack of normal stars with [Fe/H] > 
+0.2. This suggestion is supported by large samples of stars in the literature. There are 22% 
normal stars with [Fe/H] > +0.2 in Valenti & Fischer (2005)'s sample at the temperature 
range of 5600-5900 K, but note that their temperatures, derived from the fitting of synthetic 
spectra with observed spectra, are 107 K higher than our values for stars in common, and 
the sample is too small for statistics. With photometry-based temperatures and a larger 
sample of stars from Nordstrom et al. (2004), this fraction is reduced to 2%, indicating a 
low probability of super-metal-rich stars in the solar neighbourhood. This is probably due 
to the fact that most metal-rich stars tend to have formed planetary systems, and thus 
it is difficult to obtain a comparison sample that consists of super-metal-rich stars not 
harboring a planet. This is consistent with the current knowledge that SWPs have higher 
mean metallicities than normal stars. 

The second feature from Fig. 1 is that six of the 16 SWPs versus two of the 20 normal 
stars have undetected lithium lines. All of the six SWPs with undetected Li lines are main- 
sequence stars while the two normal stars with depleted Li abundances are subgiants, as 
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shown below (Fig. 5). It seems that main-sequence SWPs in the temperature range of 5600- 
5900 K destroy their Li much more easily than do normal stars. A Kolmogorov-Smirnov test 
is applied to stars with and without planets, and the cumulative distribution functions are 
shown in Fig. 2. It shows that the maximum value of the absolute difference between the 
two distribution functions is 0.54, and the probability for the two samples being the same 
group is 0.007. This is also evident from Israelian et al. (2004), as shown in Fig. 3, in which 
eight of 16 main-sequence SWPs at 5600 < T cS < 5850 K only have upper limits for Li 
abundances, while three of the 14 normal main-sequence stars from Chen et al. (2001) are 
depleted in Li. Note that the two comparison samples, from this work and Chen et al. (2001), 
are independently selected with six common stars. Thus, this property could be intrinsic, 
but further investigations with large samples of both SWPs and normal stars are needed. 

If this suggestion that main-sequence SWPs have a higher probability of depleting their 
lithium abundances is true, it will provide new information on stellar evolution for SWPs. We 
suspect that the stellar evolution history between SWPs and normal stars could be somewhat 
different due to the presence of a planet, which appears to affect only Li abundances, with no 
other elements showing any difference. As suggested by King et al. (1997), the presence of 
planets or associated circumstellar disks may affect a parent star's initial angular momentum 
and/or subsequent evolution. Specifically, the conservation of angular momentum in the 
protoplanetary disk may induce increased mixing by causing rotational breaking in the host 
stars during the pre-main-sequence evolution. This effect, once it happens, will destroy the 
Li completely. Considering this suggestion as one possibility, Israelian et al. (2004) proposed 
a second hypothesis: that migration triggers tidal forces and creates a shear instability that 
leads to a strong depletion of Li abundance. Both mechanisms occur during the pre-main- 
sequence evolution, and thus they can be used to explain our result. Further theoretical 
and observational studies are very desirable to investigate whether the two mechanisms are 
reasonable. 



3.2. Li Aabundances versus T e g 

Fig. 4 shows Li abundances versus temperatures for all stars. It seems that Li abundance 
generally increases with increasing T c g, and there is no significant gap between SWPs and 
normal stars at a given temperature. Fig. 5 shows the stellar positions of the two samples 
of stars in the HR diagram as compared with the evolutionary tracks of VandenBerg et al. 
(2000) with [Fe/H] ~ +0.12 dex. Stellar ages are also estimated from this comparison, and 
they are presented in Table 1. The error of the ages due to the uncertainties of T ef[ , M v , 
and [Fe/H] is about 15%. It is clear that our two samples of stars are generally located in 
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the same regions of the HR diagram, and during the selection all these stars are classified as 
main-sequence stars with similar parameters. However, the two stars from the comparison 
sample showing depleted Li abundances, as indicated by downward-directed arrows in Fig. 5, 
are actually subgiants. Other possible subgiant stars are not evolved far away from the main 
sequence, and there is no difference in Li abundance between SWPs and normal stars. 

This observation that subgiants tend to have depleted Li abundances reminds us to 
investigate the stellar positions of SWPs at 5600 < T c g < 5850 K in Israelian et al. (2004). 
Here we redetermined the temperatures for SWPs with 5600 < T e ^ < 5850 K from Israelian 
et al. (2004) in the same way as Chen et al. (2001). In Fig. 3, the stellar positions of these 
SWPs with new temperatures are plotted against normal stars from Chen et al. (2001) in 
this temperature range with [Fe/H] > —0.3. It is clear that a few SWPs in Israelian et 
al.'s (2004) sample based on the updated temperatures are actually evolved subgiants, and 
downward-directed arrows in Fig. 3 indicate stars with undetected Li I A6708 lines. The fact 
that evolved subgiants tend to have depleted Li abundances is consistent with the knowledge 
of stellar evolution theory. 

It has been suggested in Ryan (2000) that the Li abundances of subgiants can be some- 
what depleted due to the dilution proces without complete destruction, as opposed to main- 
sequence stars. In this scenario, once a star evolves onto the subgiant branch, it dilutes 
surface Li when convection brings Li-poor material from deep layers to the surface. This 
dilution is not obvious for slightly evolved subgiants in the present study, and a few subgiants 
with 1.2 M show higher Li abundances than main-sequence stars. Moreover, there is no 
difference in Li abundance between planet hosts and normal stars on the subgiant branch. 
Actually, the dilution process mentioned above is difficult to detect, and its effect is masked 
by the large scatter of Li abundances for even slightly evolved subgiants. It is impossible to 
investigate how the presence of a planet will affect lithium evolution for evolved subgiants 
in the present work, since their Li abundances have already been depleted due to their low 
temperatures. 

As shown in Fig. 4, there is a substantial scatter of Li abundances at T eff ~ 5800 K. 
Considering the uncertainty in temperature, it is difficult to investigate whether the scatter 
is real. There seems to be a tendency that stars in the upper envelope of the Li versus 
T e ff diagram seem to be slightly younger than stars with the lowest Li abundances in the 
T e ff < 5800 K range. However, since there is a lack of young stars in the low-temperature 
range and of old stars at the high-temperature edge in our sample, we refrain from drawing 
any firm conclusions. A larger sample of stars with different ages at a given temperature is 
needed to clarify this issue. 
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3.3. Li Anomalies of Stars With Planets and the Planet Formation 

The excess of Li depletion in SWPs in the temperature range of 5600 < T c g < 5850 K 
suggested by Israelian et al. (2004) is based on the fact that Li abundances for SWPs are lower 
than for comparison stars. A closer comparison shows that the average of Li abundances 
for 24 SWPs at this temperature range in Israelian et al. (2004) is < 0.99 ± 0.45, while this 
value is 2.08 ±0.49 for the comparison stars from Chen et al. (2001). In our work the average 
Li abundance is < 1.65 ± 0.57 for SWPs versus 2.15 ± 0.49 for the comparison stars. As 
presented above, this anomaly in Li abundance of SWPs might reflect the effect of planet 
presence on stellar evolution and cannot be used to provide constraints on planet formation. 
In addition, this effect is difficult to detect from stars in the subgiant branch or later stages 
if their Li abundances have already been depleted by other mechanisms. Moreover, the 
deviations of the average Li abundance between stars with and without planets are 0.27 dex 
for T cff > 5800 K and 0.58 dex for T cS < 5800 K. It seems that Li anomalies concentrate in 
the range of T eff < 5800 K, and thus they do not relate to the Li-dip in the Hyades, which 
happens in the temperature range of 6300-6600 K as first found by Boesgaard & Tripicco 
(1986). 

In connection with planet formation mechanism, our result does not support the ac- 
cretion scenario of planet formation. Theoretically, Montallban & Rebolo (2002) predicted 
quantitative estimates of the main-sequence evolution of stellar surface lithium after planet 
ingestion. They found that the preservation of 7 Li occurs in a large mass range of 0.9-1.3 
M & at solar metallicity. For SWPs with a typically super-solar metallicity, the mass range is 
slightly reduced. According to their theory, if a 1 Mj planet were dissolved in the convective 
envelope of an SWP at a sufficiently old age for extra-mixing to be inefficient, it would pro- 
duce an enhancement of 0.3 dex in Li abundance, which is unfortunately within the scatter 
of Li abundance at the temperature range we investigated. But Li abundances for stars with 
T e ff > 5900 K in Israelian et al. (2004) do not show such a large scatter, and the similar 
abundance between SWPs and the comparison stars argues against the accretion scenario. 
Furthermore, most SWPs known so far have giant planets with masses larger than 5Mj, 
and if 10Mj were engulfed the enhancement would reach 0.9 dex, which is easily detected. 
Therefore, we favor the scenario that the surface composition of SWPs is not polluted by 
the accretion of planetary material, and the high mean metallicity of SWPs is primordial. 

4. Conclusions 

We have made a parallel abundance analysis for 16 SWPs and 20 normal stars at the 
metallicity of -0.2 < [Fe/H] < +0.4 and the temperature range of 5600 < T cff < 5900 K. 
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Our results show that there is a higher probability for SWPs to have depleted Li abundances 
than for comparison stars, which confirms the possible Li abundance anomalies for SWPs at 
this temperature range reported by Israelian et al. (2004) after excluding evolved subgiants. 
We proposed that the presence of planets around host stars may affect stellar evolution by 
inducing additional mixing or shear instability, which leads to a high probability of destroying 
surface lithium while other elements are not affected. This effect can only be detected for 
main-sequence SWPs since other depletion mechanisms would not have started to act on 
main-sequence stars. In view of this, the effect of planet presence should be considered 
when the large scatters of the Galactic lithium in field disk stars at the temperature range 
of 5600 < T cS < 5900 K are investigated. 
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